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We present high resolution two-photon excitation microscopy studies 
combining two-photon fluorescence (TPF) and second harmonic generation 
(SHG) in order to probe orientational distributions of molecular ensembles 
at room temperature. A detailed polarization analysis of TPF and SHG 
signals is used in order to unravel the parameters of the molecular orienta- 
tional statistical distribution, using a technique which can be extended and 
generalized to a broad variety of molecular arrangements. A polymer film 
containing molecules active for TPF and/or SHG emission is studied as a 
model system. Polarized TPF is shown to provide information on specific 
properties pertaining to incoherent emission in molecular media, such as 
excitation transfer. SHG, being highly sensitive to a slight departure from 
centrosymmetry such as induced by an external electric field in the medium, 
complements TPF. The response of each signal to a variable excitation pola- 
rization allows investigation of molecular behavior in complex environments 
which affect their orientations and interactions. 

I. Introduction 

Multi-photon microscopy has amply demonstrated its assets in the study of 
broad variety of physical and biological phenomena. Among the effects aris- 
ing from two-photon excitation, two-photon fluorescence (TPF) has been at 
the focus of much attention for many reasons, fl] Contrary to linear fluores- 
cence, for which background noise rejection requires the implementation of 
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confocal detection, TPF exhibits an intrinsic high spatial resolution due to 
the built-in quadratic power dependence on the excitation intensity. In addi- 
tion, the use of infrared excitation minimizes both optical damage and back- 
ground scattering in complex samples such as living cells. ;2 This method has 
shown significant advantages in imaging of complex biological environments 
down to the single molecule level. [3] In parallel, developments in molecu- 
lar engineering have lead to the design of chemical probes exhibiting very 
high two-photon absorption cross sections. [U |5] Such results have triggered 
important advances in the field of two-photon imaging and other processes 
such as optical limiting. Finally, multi-photon microscopy techniques have 
been extended to a larger variety of nonlinear optical processes such as 
three-photon induced fluorescence, as well as second- and third-harmonic 
frequency generation both in far-field and near-field microscopies. [HI E[|HJ|!i] 

Second harmonic generation (SHG) is the result of coherent emission 
from anharmonic oscillators upon two-photon excitation. At macroscopic 
scale, SHG emission vanishes in centrosymmetric media, since the x re_ 
lated nonlinear susceptibility is an odd-rank tensor. In molecular media, 
non-centrosymmetric ordering is traditionally obtained by electric field pol- 
ing of intrinsically non-centrosymmetric molecules, such as dipolar 7r-conju- 
gated systems functionnalized with adequate acceptor and donor groups. ^U] 
As SHG is highly sensitive to even a slight departure from centrosymmetry in 
molecular media, this effect has been widely used to study surface and inter- 
face properties in physics and chemistry, with monolayer sensitivity. [Ill I12j 
SHG signals have also been measured in artificial vesicles. [131 live cells. [14] 
and in biological membranes in the presence of an electric potential. |15j 

In this work, we study TPF and SHG in model molecular systems and 
take advantage of the specific properties and complementarities of these 
nonlinear processes. The main difference between them lies in the coherent 
nature of the emitted signals and their relative spectral features. Contrary 
to fluorescence, which is based on the time-delayed relaxation of a radiative 
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level, off-resonant SHG occurs instantaneously upon excitation. Further- 
more, the fluorescence emission is affected by a Stokes shift, whereas SHG 
occurs at half the incident wavelength and is phase-correlated with the exci- 
tation field at the fundamental frequency. These specific features allow one 
to distinguish them either spectrally or temporally. |16| 

The main advantage of combining the two processes lies in their respec- 
tive anisotropy sensitivity to molecular orientation under polarized excita- 
tion. In a centrosymmetric medium, SHG vanishes while TPF is always 
present and can be anisotropic if the molecular orientational distribution 
symmetry is axial. [T7] Moreover, SHG appears in a polar medium, and is 
intrinsically anisotropic. The difference between these two effects is due to 
the distinct symmetry properties of the tensorial susceptibilities that are 
involved, TPF and SHG being described by even and odd order tensors re- 
spectively. The related macroscopic susceptibilities reflect both molecular- 
scale and macroscopic-scale arrangement properties, such as orientational 
distribution. Measuring TPF and SHG anisotropics is therefore a direct 
way to probe a given molecular organization in a molecular ensemble that 
is either naturally ordered, or that has been subject to the application of 
an external perturbation, such as an electric field. 151 119j a coherent com- 
bination of optical excitations, 20 or a combination of these. |2*T] Complex 
tensorial distributions can be unravelled using such approaches. |22| 

In order to fulfill the symmetry requirements discussed above, we chose 
a model medium consisting of an amorphous polymer matrix doped with 
fluorescent and/or dipolar nonlinear active chromophores. Macroscopic cen- 
trosymmetry breaking of the medium is induced by electric field orientation 
of these chromophores, using planar geometry suited to the two-photon mi- 
croscopy experimental setup. Applying an external electric field results in 
a controllable external degree of freedom used to monitor the molecular 
order symmetry. Furthermore, the reflection geometry of the inverted mi- 
croscopy configuration used in this work allows one to measure both TPF 
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and SHG signals without any dependence on the coherence length of the 
nonlinear field propagation. The experimental configuration presented here 
can be easily applied to a variety of situations such as in biological mem- 
brane studies where electric potentials might induce the orientation of a 
molecular probe. Such studies can be potentially down-scaled to a low num- 
ber of molecules and even to single macromolecules, where interactions at 
the molecular level or global conformational changes can affect the collective 
emission anisotropy properties. 

We present a detailed analysis of the orientational molecular ordering in 
the sample using both TPF and SHG processes. We show that the rotating- 
polarization analysis can bring additional informations that might be hid- 
den in ratiometric anisotropy measurements. This analysis can be straight- 
forwardly adapted to a large variety of environments, and to dynamical 
effects. [TH] Accounting for the instrumental features of inverted microscopy, 
with high numerical aperture (N.A.) objectives collecting light at wide angles 
of emission, we show that a simple model allows one to retrieve information 
on localized dipole orientations. As we demonstrate in this work, the ana- 
lysis of nonlinear anisotropies requires therefore particular attention to the 
polarization state at the focus of a high N.A. objective and the definition of a 
calibration procedure. In the following sections, we discuss the optical prop- 
erties of both TPF and SHG processes in random and ordered media, taking 
into account the influence of excitation transfer between chromophores. 

II. Experimental section 
A. Sample Preparation 

The samples, studied at room temperature, contain active chromophores, ei- 
ther 4-dicyanomethylene-2-methyl-6-(p-(dimethylamino)styryl)-4H-pyran (DCM) 
or 4-(N-ethyl-N-(2-hydroxyethyl))amino-4'-nitro-azobenzene (Disperse Red 1, 
DR1). The DCM chromophore exhibits a high fluorescence quantum yield 
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and a non negligible molecular quadratic polarizability (3. 23 The DR1 chro- 
mophore, which has been extensively studied for its intrinsic high quadratic 
nonlinear efficiency. |24| is not fluorescent because of quenching by fast pho- 
toisomerization. Both chromophores exhibit a large permanent dipole mo- 
ment which enables electric field orientation. 

The molecules are dispersed with a maximum concentration of 10% 
by mass in a high T g (120°C) polymer matrix of polymethyl-methacrylate 
(PMMA), which is not cured after deposition. The polymer films are then 
deposited by spincoating either on a microscope slide or on an electrode 
substrate. |25| In order to achieve electric field alignment in the plane of the 
sample, we use a set of gold transverse planar electrodes of 50 nm thickness, 
separated by a gap of 10 /xm, as represented in Figure ^ The electrodes 
are fabricated by photo-lithographic patterning of UV sensitive photoresists 
on the gold layer, which is deposited by sputtering on a microscope slide 
covered with a 100 nm PMMA adhesion layer. The thickness of the active 
PMMA-dye layer is of the order of 100 nm, in which the active molecules are 
expected to orientate perpendicular to the electrodes in the sample plane. 

B. Experimental Setup 

The experimental setup used for nonlinear microscopy is described in Fi- 
gure ^ The source for nonlinear excitation is a mode-locked Ti:Sa laser 
which produces 120 fs pulses at a fundamental wavelength of 987 nm with a 
80 MHz repetition rate. The laser beam is focused on the sample by a high 
numerical aperture oil immersion objective, leading to a spatial resolution 
of 400 nm. Typical incident energies range from 0.01 nJ to 0.1 nJ per pulse. 
The TPF and SHG signals arising from the sample are collected by the same 
objective, and then directed to a polarizing beamsplitter and a set of two 
avalanche photodiodes operating in the photon counting regime. We select 
either the SHG signal or the TPF signal by appropriate interference filters. 
The spectral distribution of the emitted light can be analyzed in parallel, 
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Figure 1. Nonlinear microscopy setup: HP filter: high pass filter which 
rejects visible light; P: polarizer; HWP: rotating half-wave plate; Obj: mi- 
croscope objective (x60, NA=1.4); DM: dichroic mirror; F: interference filter 
selecting either the SHG signal at 493 nm or the TPF signal at 580 nm; IR 
filter: filter rejecting the residual incident (near IR) laser beam; PBS: po- 
larizing beamsplitter selecting X and Y polarization states of the emission; 
APD: avalanche photodiodes. An additional beamsplitter (BS) can be in- 
troduced for spectroscopic measurements. During the polarization analysis, 
we remove this beamsplitter in order to not perturb the polarization state 
of the detected signal. 
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using a spectrograph coupled to a highly sensitive CCD camera. 

For both TPF and SHG, the polarization analysis consists of rotating 
the incident polarization from 0° to 360° and recording the corresponding 
emissions on the two perpendicular analyzed polarization directions. The 
geometry of the system is schematically represented on Figure ^ The Z 
axis is along the optical axis, perpendicular to the sample. The X and Y 
axes, lying in the sample plane, provide an analysis framework defining the 
polarization directions detected by the two photodiodes. They also coincide 
with the p and s incident polarizations on the dichroic mirror. The incident 
excitation polarization at the fundamental frequency lo can be written at 
the focus point as:Pfi] 



E 

E(a, 5, 7, cot) 



cos a cos(tot) 
(1 — 7) sin q cos(o;i + 5) 



where E is the field amplitude and the rotating angle parameter a defines 
the incident polarization angle in the (X,Y) framework. In this model, we 
introduce two external parameters, 7 and 5, to account for polarization mix- 
ing effects in the excitation/detection setup. The parameter 7 represents the 
difference of p and s reflectivity from the dichroic mirror. The phase shift 
5 represents the ellipticity appearing after reflection on the dichroic mirror 
with an incidence angle of 45°. The X(p) and Y(s) polarizations have been 
observed to be unaffected by any ellipticity, allowing us to assume that ellip- 
ticity is only introduced for intermediate a polarization angles. Note that 
in this model, we neglect the longitudinal Gouy phase shift factor Z contri- 
bution of the incident field. |27| the sample thickness being much lower than 
the coherence length of the nonlinear process. |281l29j After passing through 
the dichroic mirror, the light transmitted (reflected) by the polarizing beam- 
splitter used for analysis is polarized in the Y (X) direction, corresponding 
to Iy (Ix) detected intensities (see Figure ^). 
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III. Results and Discussion 

A. Two-Photon Fluorescence Microscopy 




Figure 2. Fluorescence intensities Ix (in dark grey) and ly (in light grey) 
emitted by an isotropic assembly of DCM molecules, shown in a polar plot 
as a function of the a angle varying from 0° to 360° . The integration time is 
100 ms per point. The four diagrams have been obtained for four different 
concentrations in mass of DCM in PMMA: (a) 10% (b) 0.83% (c) 0.085% (d) 
0.022%. The continuous lines represent the fit according to the expressions 
given in Appendix A. The four graphs have been normalized to the same 
intensity. The distinction between the two intensity amplitudes Ix and ly in 
(a) is due to a slight difference in detection efficiency in the two polarization 
channels 
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In order to calibrate the polarization response of the nonlinear microscopy 
setup and to account for the effect of the high numerical aperture objec- 
tive on the light polarization, we first study the fluorescence emitted by an 
isotropic molecular distribution. The sample used for this purpose is a guest- 
host polymer matrix of PMMA embedded with the fluorescent dye DCM at 
different concentrations. Figure [21 shows the experimental variations of Ix 
and Jy fluorescence intensities with respect to the polarization direction a 
of the excitation beam. In order to explain the observed features, several 
experimental considerations have to be accounted for, in addition to the 
incident field polarization parameters detailed previously. 

First, the high numerical aperture of the objective affects the polari- 
zation radiated by the dye molecules. [2*7] This effect can be modelled us- 
ing a calculation similar to Ref. |3U| relative to one-photon fluorescence. 
Indeed, TPF emission is still occurring from a one-photon allowed transi- 
tion, independently of the excitation pathway. We define the direction of 
the emission dipole of the rod-like DCM molecule by a set of two angles 
(9,(f)) = The fluorescence intensities Jx and Jy emitted by this single 
dipole p(Q) = (nx(Q), fJ>Y (^)j Hz(&)) and measured by the two detectors 
can be written as expressed in Appendix A: 

Jx(fi) = k x i4c(n) + k 2 vUn) + k 3 

Jy (O) = K 2 fx («) + My («) + K 3 ix\ (O) (2) 

where 7^1=2.945, 7^2=0.069 and K 3 =1.016 in the case of a NA=1.4 objective 
such as the one used in the present work. K\, K 2 and K 3 represent the 
mixing of polarization components in the emission, as a consequence of the 
collection of light at very wide angle. 

Second, the two-photon nature of the excitation process has to be taken 
into account. Contrary to a 1-photon excitation which depends on the square 
of the incident field E, the 2-photon excitation probability of a dipole /2 is 
proportional to |/2./5| 4 . Therefore, the polarization state of the incident IR 
beam plays a crucial role. We will furthermore suppose that the emission 
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and excitation dipoles are parallel for DCM molecules. |31j 

Third, we assume a normalized molecular orientational distribution func- 
tion /(f2) which is set equal to 1/4-7T in the case of an isotropic molecular 
ensemble in the polymer matrix. Using the previous assumptions, the re- 
spective fluorescence intensities in the X and Y analysis directions can be 
expressed as: 

^(direct) = J m j\4j l{n )f(n)dn (3) 

after orientational averaging with df2 = sin 9 d8 dcj) and subsequent temporal 
average represented by the (...) notation. 

Finally, Figure EJ shows that the polarization response depends strongly 
on the molecular concentration, which is the signature of a possible excita- 
tion transfer between fluorescent molecules. The effect of such transfer in 
highly concentrated media is to decouple the excitation step from the emis- 
sion one, j32 resulting in a depolarization process, which leads to the same 
signal in the two analyzing channels. This is indeed what is observed in 
Figure The excitation transfer efficiency between neighboring-molecules 
depends on the relative orientation factor k 2 between two neighboring chro- 
mophores, the "donor" p%(9i, <j>\)) and the "acceptor" /U 2 (^2! fa, p, £, #i, </>i)):_33 

K 2 {p, @2, 4>2) = (2 cos p. sin 02- cos 4>2 + sin/3, sin 2 - sin^) 2 (4) 

In this equation, (9%, (j>i)=Qx defines the orientation of the donor pi, (p, £)=0 
defines the orientation of the vector ui 2 connecting the two chromophores in 
the (#1, 4>\) framework and (62, 4 > 2)=^2 defines the orientation of the accep- 
tor p2 in the (p, £) framework. The expression of the vector /2 2 (0 2 , (f>2,P, £> #1, 4>i) 
in the macroscopic framework (X, Y, Z) is given in Appendix B. The fluo- 
rescence intensity after energy transfer can then be written as: 

^TPF Ctransfer) = f f f \^ (^.E^K 2 (p, 2 , <£ 2 ) Jl(?h, fi X , 0)/(0 1 )/(0 2 )dfidOldfi 2 
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where Ji{Q,2 1 fii, is given in Appendix B. The total fluorescence intensity 
can then be written as: 

jTPF _ jTPF (direct) _j_ j- jTPF (transfer) /g\ 

where T is a fitting parameter which quantifies the transfer rate. As seen 
in Figure |^1 this model is in good agreement with the data, using the ex- 
perimental factors 5=1.15 rad for the field phase shift responsible for the 
incident ellipticity and 7=0.02 for the dichroism. These parameters have 
been measured separately by ellipsometry at our fundamental wavelength. 
The effect of ellipticity is to give to the polarization response a cross shape in- 
stead of the expected uniform circle in the case of randomly oriented isolated 
molecules. Figure corresponds to a high concentration sample resulting 
in a complete depolarization of the emission for which a transfer rate T=l 
allows one to fit the data. Due to the intermolecular long-range coupling 
between the two dipoles, the energy transfer rate T depends mainly on the 
inter-molecular distance r in the polymer and is proportional to 1 /r 6 . r was 
varied by decreasing the DCM concentration in mass, namely from 10% to 
0.83% and 0.085%, down to 0.022%. In each case, an average intermolecular 
distance can be estimated, giving values of 6.0 nm, 7.0 nm, 8.9 nm and 10.9 
nm respectively. The effect of a decreasing concentration appears clearly 
in Figure |5J with the lateral lobes of Ix and Iy disappearing progressively 
in the Y and X directions respectively. Note that for the concentration 
corresponding to Figure the excitation transfer is negligible, and each 
detector detects mainly its preferential polarization direction as expected 
for an incoherent optical process. The energy transfer rates obtained from 
the fits of Figure|21are respectively 1, 0.15, 0.03 and 0.02. These values are 
represented on Figure The dependance with the intermolecular distance 
r is compared with a 1/r 6 power law, showing a qualitative good agreement. 
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Figure 3. Log-log plot of the experimental transfer rate versus estimated 
intermolecular distance, as the DCM concentration in the polymer matrix 
is varied. A 1/r 6 power dependance is also plotted as a guideline. 



B. Second Harmonic Generation Microscopy 

The same setup can also be used to study the influence of a non-centro- 
symmetric contribution of the molecular distribution on the optical response 
upon two-photon excitation. In this case, a SHG signal arises from the po- 
lar orientation of nonlinear active molecules. For this purpose, the DCM 
doped polymer matrix is spincoated on a set of two transverse electrodes, 
as described in the first part of this paper. A high voltage of about 300 
V, corresponding to a poling field of 3 x 10 7 V.m , is applied by means 
of electrical contacts on the electrodes. As the DCM molecules have a non 
negligible permanent dipole moment of about 10 Debve. [23"] they will have 
a tendency to lean towards the electric field direction, thus creating a non- 
centrosymmetric distribution. Moreover, due to the nonlinear susceptibility 
j3 of DCM, [22! the molecular orientation can be directly monitored by the 
onset of a SHG signal in addition to fluorescence as reported in the previous 
section. In the present configuration, the SHG signal is detected in a reflec- 
tion configuration. Since the thickness of the molecular layer is much smaller 
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Figure 4. (a) Two-photon emission from a polymer film doped with DCM. 
A SHG peak appears next to the fluorescence band when a high voltage is 
applied to the sample (dashed line: V; continuous line: 300 V). (b) SHG 
modulation from DCM molecules under an ON-OFF poling field. The high 
background is due to the DCM residual TPF, which lies within the filter 
bandpass. Note that it would be possible to avoid such a signal by working 
off resonance for the SHG. (c) Two-photon emission from a polymer film 
doped with DR1, showing only SHG. The width of the SHG peak is limited 
by the incoming IR femtosecond laser linewidth. (d) SHG modulation for 
DR1 molecule. In both cases, the incident IR laser beam energy is about 
0.1 nJ per pulse. 



13 



than the coherence length for phase-matching, the molecules can be simply 
considered as a coherent ensemble of nonlinear radiating dipoles. The SHG 
signal therefore appears with the same intensity either in the transmitted 
or in the reflected direction. As shown on Figure 0k, the overall two-photon 
emission spectrum exhibits a narrow SHG peak at 2u>, emerging in the tail 
of the broadband fluorescence when the electric field is applied. The SHG 
and TPF signals can be easily selected by introducing appropriate optical 
filters. Figure shows the temporal evolution of SHG upon the ON-OFF 
application of the electric field: since the SHG signal disappears when the 
electric field is turned off, on can assume that the molecules still have a 
high degree of mobility in the polymer matrix in spite of its high T g . This 
mobility is probably a result of the broad size distribution of cavities in the 
polymer with a tail in the large size limit. 

When spectrally selecting the TPF signal only, the polarization response 
of the sample does not undergo significant changes in the presence of the 
static electric field. The fluorescence sensitivity to molecular orientation can 
be explored by calculating similar polarization patterns as those observed 
in Figure [21 accounting for the molecular orientational distribution induced 
by electric field orientation. The calculated fluorescence patterns, shown on 
Figure [51 have been computed considering a permanent poling field Eq and 
using the Boltzmann distribution function: 

(jl(n).E Q )/k B T) 

/ rn) = — = (7) 

j e (fi(n).E )/k B T) d Q 

where ks is the Boltzmann's constant and T the poling temperature. As 
seen from Figure 03 in our experimental conditions, i.e. with an electric 
field of 3 x 10 7 V.m -1 and at room temperature, the fluorescence patterns 
exhibit only slight differences compared to an isotropic molecular distribu- 
tion. This result is not surprising since the TPF intensity is proportional to 
the number N of active molecules per unit volume, whereas SHG intensity 
is proportional to N 2 as a result of the coherence of this radiating process. 
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Moreover, the weak effect of the electric field on the fluorescence anisotropy 
is also due to the broad statistical orientational Boltzmann distribution at 
room temperature. Figure shows that orientational effects should be 
strongly increased with higher electric fields, which are however difficult to 
reach due to electric breakdown. 




Figure 5. Influence of the poling field Eq, aligned along the X axis, on 
polar polarization plots of the calculated TPF intensities Ix (dashed line) 
and Jy (continuous line) emitted by an assembly of fluorescent molecules. 
Energy transfer effects are not taken into account in these calculations. All 
graphs have the same scale (a) Fluorescence diagram for an isotropic mole- 
cular distribution, (b) Fluorescence diagram with a field of 3 x 10 7 V.m -1 , 
such as used in our experiment. There is no clear difference with the first di- 
agram, (c) The effect of the applied field becomes visible for 3 x 10 8 V.m -1 , 
corresponding to a voltage 10 times higher than the experimentally used 
values. 

The SHG signal is therefore a better suited probe of the field-induced 
anisotropic molecular distribution. In order to retrieve SHG information 
independently from fluorescence properties, we choose to study more specif- 
ically samples doped with the DR1 molecules, so as to prevent from possibly 
misleading combinations of SHG and TPF signals. This situation is akin to 
using DCM molecules with non-resonant excitation, where TPF would then 
be inefficient. DR1 has a non-negligible dipole moment of 8.7 Debve. [21] 
and a molecular hyperpolarizability /?=289xl0 -40 m 4 V _1 at zero frequency, 
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150 counts 800 counts 

Figure 6. Polar plots showing the SHG signals (integration time: 100 ms) 
Ix (in dark grey) and ly (in light grey) emitted by an assembly of DR1 
molecules poled under an electric field of 3 x 10 7 V.m -1 as a function of the 
incident polarization angle varying from 0° to 360°. The different diagrams 
have been obtained for different orientations of the electric field, as indicated 
on the Figures, (a) and (b): Ix and ly for Eq along X, (c) and (d): Ix and 
ly for Eq along Y, (e): Ix and Iy for Eq at 45° of X and Y. The continuous 
lines represent the fit according to the expressions given in Appendix C. 
The situations where Eq is along the X or Y axis are symmetric, whereas a 
rotation of 45° mixes the nonlinear polarization responses together. 

which amounts to twice that of DCM.j52] Figures and^l show the ef- 
fect of electric field poling of DR1 using the same electrode system . The 
same polarization analysis is applied to investigate the SHG responses, and 
Figure shows the evolution of the SHG signal with respect to the rotation 
of the incident polarization a. 

In order to fit the data, we developed a model elaborating on Ref. 
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^IIM], based on the macroscopic second order susceptibility x (2w; uj, uj), 
which is related to the molecular hyperpolarizability P(2uj;uj,uj) tensor and 
the molecular distribution function /(O). Earlier models are extended herein 
so as to account for specific features pertaining to the two-photon microscopy 
setup, as detailed in Appendix C. The rod-like DR1 molecule is composed 
of a conjugated tt system with an electron acceptor group at one end and a 
donor group at the other end. The nonlinear susceptibility (3 has therefore 
only one non-zero component (3 ZZZ , where z defines the molecular axis direc- 
tion. We denote £1 = {9, 4>) the Euler angles defining z in the macroscopic 
framework. With N molecules per unit volume, the susceptibility tensor 
components can be expressed as: 

Xuk = N(3 ZZZ J cos(z, /) cos(z, J) cos(z, 2f)/(fi)dfi (8) 

where the indices I,J,K = X,Y,Z are the coordinates in the macroscopic 
framework and f(£l) is the molecular orientational distribution function 
given by Eq. (J7|). The cos(z, I) functions are the (f2) angle dependent pro- 
jections of the z axis on the / axis. For symmetry reasons, a non-zero value 
for x requires a macroscopic polar order, which is imposed here by the 
static electric field Eq applied in the plane of the sample. From the Xijk 
tensor coefficients, we can infer the induced nonlinear polarization p( 2 )(2w) 
defined as: 

PP{2u) =Y,X?LEj(")Ek(cj) (9) 

J,K 
(2) 

Note that this expression of P) (2u>) does not take into account the local 
field factors ff^, fj and /^-, which is reasonable in the present context of 
relatively low molecular concentration and weak poling strength. As de- 
tailed in Appendix C, this polarization allows us to compute the 2uj radiated 
field. The correction factors accounting for the collection of light at wide 
angles are estimated with a model similar to the one used for TPF. Since 
SHG is a coherent process arising from the induced dipole radiation, the 
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collected field can be written as: 

£? hg {2uj) = b IXX E\ + b IYY E Y + bixyExEy (10) 

where bjjn are constant coefficients defined in Appendix C. The resulting 
intensity is proportional to \£ (2u>)\ 2 , and the SHG intensities /| HG de- 
tected by the photodiodes after temporal averaging can be finally expressed 
as: 

lf^ Y = bj xx 1% + bj YY Ey~+(2 bixx birr + bj XY ) E\E\ 

+ 2 b IX x bixY E\E Y + 2 biw bixv E X E\ (11) 

Figure El represents the experimental data and the theoretical fits for several 
orientations of the poling field Eq relative to the (X, Y) axes, showing a 
good agreement. These fits also take into account the ellipticity parameter 
5 = 1.15 rad and the 7 factor of 0.02 determined previously. Contrary to 
fluorescence, the low intermolecular distance in the polymer matrix is seen 
to not affect the polarization response of SHG. These results show that a 
thorough SHG polarization analysis is able to account for a slight change 
in molecular organization. In particular, the symmetry features of SHG po- 
larization polar plots shown in Figure are consistent with the rotation of 
the electric field in the (X, Y) plane. Such polarization signatures can be 
exploited to identify specific molecular orientation directions in the sample 
plane. In the present case, the 2oj optical response is spatially uniform over 
a raster scan of 5 /xmx5 /im (data not shown). This provides clear evidence 
of the good homogeneity of the field orientation between the two planar elec- 
trodes. Such studies are being currently extended to the mapping of mole- 
cular polar-orientation in electro-optic devices based on complex electrodes 
designs, or to lower-scale changes in molecular orientation patterns such 
as those observed in Langmuir-Blodgett films. [37] 
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IV. Conclusion 



Polarization analysis of the emission response from an assembly of molecules 
under two-photon excitation contains significant information on their orien- 
tational distribution. Focusing on either TPF or SHG furthermore allows 
one to study either the incoherent, axial even-order processes, or the coher- 
ent polar odd-order processes. Combination of these two techniques provides 
a complete set of informations regarding possible excitation transfer between 
molecules or orientation of a collection of molecules with very accute sensi- 
tivity to the onset of non-centrosymmetric orientational patterns. 

The approach presented in this paper can be applied to a broad vari- 
ety of molecular media, where unknown orientational distribution functions 
are expressed using spherical harmonics which can be singled-out via their 
symmetry properties. By taking advantage of the high optical resolu- 
tion provided by two-photon microscopy, interesting perspectives can be 
expected in the study of structures involving sub-microscopic scale effects 
such as in nanocrystals 3]3 0H1 and nanostructured media supporting local 
field enhancements. [lU] 

Appendix A: Polarization response of TPF in a high NA mi- 
croscope setup for wide angle fluorescence collection 

The following model is developed for comparison between polarization re- 
sponses in two perpendicular polarization directions (see Figure In the 
subsequent derivations, we omit all efficiency parameters that may appear 
in the intensity expressions, since such parameters are the same for each 
channel of polarization direction. 

In order to express the intensity of fluorescence emitted by a single dipole 
set at the focal point of the microscope objective, we consider the radiation 
diagram emitted by this dipole /2, which orientation is defined by the angles 
(6,4>), as indicated on Figure [Tj The far field radiated in the direction of 
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z 




Y 



Figure 7. Notations introduced to account for the wide angle collection of 
the emitted light by a dipole, located at the origin of the (X, Y, Z) frame- 
work. The direction of the dipole jl is defined by the angles (0,(j>). The 
direction of the wave vector k corresponding to the radiated field, is defined 
by the angles (u, v). 



the wave vector k(u, v) is: 

£*TPF (radiated) ^ A (jl A k) (12) 

were the proportionality factor contains constant radiation field terms that 
are not used in the following calculations. As the orientation of the wave 
vector k is only defined by the angles (u, v), the vector £ TPF ( radiated ) can 
be expressed as: 

£-TPF (radiated) ^ ^ ^ = ^ ^ ^ jj^ y)+ ^ ^ ^ ^ „) + ^( 0> „) 

(13) 

where £/i, C^, C/3 are unit vectors depending only on the parameters (u,v). 
The field transmitted by the objective can then be expressed as: 

c*TPF f7?l C* TPF (radiated) (i^ 
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where [7Z] represents the rotation matrix simulating the infinity-corrected 
objective. [R] is therefore the product of three successive rotations 3Qj (ro- 
tation of — v around Z, rotation of —u around Y and rotation of v around 
Z) so as to convert any input incidence on the objective into an output ray 
parallel to the optical axis Z, namely: 



cos u cos 2 v + sin 2 v cos v sinv(cosu — 1) — sin n cos w 
cos v sini^cosit — 1) cos u sin 2 v + cos 2 v — sin -u sin?; 
smu cosv sinti sinw cos it 



(15) 



The vector £ TPF (it, v, 9, </>) of the trasmitted field can then be expressed as: 

£x FF (u,v,9,(j)) = f x {u,v)nx{0,(j>) + fy(u,v) Hy(0,<P) + fz(u,v) nz{0,<t>) 
£y PF (u,v,9,(f)) = g x {u,v)nx(0,4>) + gy{u,v) Hy(9,(/>) + 9z{u,v) fi z (9, 0)(16) 

where fx, fy, fz, 9x, 9y and gz are functions of the (u,v) parameters. 

Since fluorescence light is emitted incoherently, the intensities detected 
by the photodetectors and coming from the single dipole f2(9,4>), are com- 
puted after integration of the square of each £ TPF component, over all the 
angles (u, v) within the half-aperture angle o bj of the objective, hence giving 
the detection probability: 



2tt 



'obj 



Ji=x,y(0,<I>) = / {£j VF (u,v,9,$)) sinududu (17) 
Jo Jo 

In the present work, the oil-immersion (n=1.5) objective has a numerical 
aperture NA=1.4. The half-aperture angle o bj is therefore equal to 1.204 
rad (69°). After integration, the previous expression reduces to: 

Jx (9, </>) = K x fi x (9A)+ K 2 &(9, <j>) + K 3 n\{9, 0) 
Jy(9,<l>)=K 2 fix(9,0) + K 1 fi Y (9, ( j>) + K 3 iJ 2 z (9,cl)) (18) 
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with: 

r-2-K r 0nhi /-27T 



K\= I I f x (u,v) sin ududv = / / (u, t>) sin ndti du = 2.945 
Jo Jo Jo Jo 

r-2n r8 ohi |-2tt rO ohi 

if 2 = / / f Y (u,v) sinudndt' = / / (u, sin udu dv = 0.069 
Jo Jo Jo Jo 



2t /'t'obj f^TT fdobj 

2/„. „.\ „:„ „, J„. J„, / / „2 



K3 = / / z (n, sinu du dw = / / ^(ti, u) sin udu df = 1.016 (19) 
Jo ' Jo Jo 

Moreover, the two photon excitation probability, which is proportional 
to |/2(f2)..E| 4 , can be easily evaluated. Since fluorescence is an incoherent 
process, we directly integrate the dipole response, product of the excita- 
tion and detection probability, over all dipole orientations. Assuming that 
the absorption and emission dipoles of each chromophore are parallel, the 
detected intensities can then be expressed as: 

iI=Iy = f \m)-E\ 4 Mtyf(n)dn (20) 



Developing |/x(S7)..E| 4 in the previous expressions, the TPF intensities, de- 
pending only on parameters (a, 5, 7), finally reduce to: 

lJ=x iY ( a , s > 7) = ^2 fiJKLM EjE K E L E M 

J,K,L,M=X,Y 



fixxxxE x + f IYYYY Ey + 6flXXYYE x E y 



+ ifixxxYE x E Y + Af IX YYYE x E Y (21) 
where the /ijklm coefficients are defined by: 

fiJKLM = J Ji(n) Hj(SI) MJf(«) l*L(to) M«) /(«) (22) 

Appendix B: Expression of the acceptor coordinates /12 in the 
macroscopic framework 

As indicated in Figure |H1 </>i) defines the orientation of the donor dipole 
pi and (p, £) defines the orientation of the vector U12 connecting the two 
chromophores in the (#i,</>i) framework. The orientation of the acceptor 
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Figure 8. Position and orientation of the acceptor fa relative to the position 
and orientation of the donor p\ . The orientation of the donor fa is defined by 
the (0\, <pi) angles. The orientation of the unitary vector uyi along the line 
connecting the two chromophores is defined by the (p, £) angles in the (0\, 4>\) 
framework. The orientation of the acceptor fa in the (p, £) framework is 
defined by the (02, $2) angles in the local framework of molecule 2 (shown 
in dotted lines). 



dipole fa is defined by the angles (#2,^2) i n the (p, £) local framework of 
molecule 2. The acceptor coordinates in the macroscopic framework depend 
consequently on the angles (02, 4>2, P, £, Q\, 4>l) an d can be written as: 



#2(02,02, P,£, 0\, 4>i) = [>l](e 1 ,0 1 )^(x,y,z)[- A ^](p if )-^(e 1 ,</ )1 ) 



(23) 



sin 02 COS (p2 
sin 02 sin $2 

COS 02 

where [A / i]( p ^)_ > (e lj ^ 1 ) represents the rotation matrix from the (p, £) frame- 
work to the (0±, (j>i) framework and [M](e 1 ,</> 1 )^(x,Y,z) is the rotation matrix 
from the (Oi,<f>i) framework to the macroscopic one. Expressions for these 
matrices are: 

sin p cos £ cos p cos £ — sin £ 
[■^](p,5H(<?i,<M = sin/9sin£ cospsin£ cos£ (24) 

cos p — sin p 
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and: 



[M\e u ^{x,Y,Z) = 



sin 9\ cos 4>\ cos 9\ cos <p\ — sin <pi 
sin 6\ sin cos #1 sin <fii cos 0i (25) 
cos 6*i — sin #1 

then be expressed as: 



The fluorescence intensity emitted by the acceptor can 

jy(o 2 ,^i,p,6 = k 2 pIx(^2, 4>2,p,£, Oi, fa) + Ki^l Y (e 2l 02, P, £,0l,<t>l) 

+ K 3 ^ 2Z (9 2 ,02,P,tOi,<Pi) (26) 
where the parameters Ki, K 2 and K 3 are defined in Appendix A. 

Appendix C: Polarization responses of SHG in a high NA 
microscope setup 

In the case of second harmonic generation, the expression of the radiated 
field is similar to the fluorescence emission, requiring mere replacement of the 
dipole moment p by the induced second order nonlinear polarization : 



£ SHG ( radiatcd )( n , v ) oc k A (P (2) A k) 



(27) 



with: 



p(2) 



r v (2) 1 

Xxxx 




" (2) " 




' (2) " 

XxXY 


y (2) 




(2) 


E Y + 2 


(2) 
Xyxy 


y (2) 




(2) 
Azyy 




(2) 
Azxy 



ExEy 



(28) 



(2) 

where the Xuk coefficients are the tensor coefficients of the macroscopic 
nonlinear susceptibility ■ The radiated field £ SHG(radiated) j g then: 



£ SHG(radiatcd) ^ , 



A-XXX A-XYY Axxy 

Ayxx Ayyy Ayxy 
A-zxx Azyy Azxy 



Ex 
El 

ExEy 



(29) 
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where the Ajjk coefficients depend only on parameters (u,v). Taking into 
account the effect of the objective (see Eq. p4[l). we can write: 





Bxxx 


Bxyy 


Bxxy 




E\ 


£ mG (u,v) = [R]£ SUG ( radiatcd ) = 


Byxx 


Byyy 


Byxy 




El 




Bzxx 


BzYY 


Bzxy 




E X E Y 



(30) 



where the Bjjk coefficients depend only on the (u, v) angles. As SHG signal 
is emitted through a coherent process, integration over all (u, v) angles in 
the objective aperture is performed before taking the square of the emitted 
field. The total collected field £ t o t HG is then obtained by integration over 
the (u, v) angular variables: 





bxxx 


bxYY 


bxxy 




E\ 


o*SHG 
°tot — 


byxx 


byyy 


byXY 




El 




bzxx 


bzYY 


bzxy 




ExEy 



where: 

/•2-7T r-e obj 

b[jK = / / Bijk(u, v) s'mududv (32) 
Jo Jo 

Since the detected intensity is proportional to |£ t ot HG | 2 ! the expression of Ix 
and Iy are therefore: 

lf=x,Y = bjxx E x + bj Y y E Y + (2 bjxx biyy + b 2 IX y) E x El 

+ 2 b IXX b IXY E\E Y + 2 b IYY bixy E X E\ (33) 
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